Exact resolution of the Baxter equation for reggeized gluon interactions 
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Abstract 

The interaction of reggeized gluons in multi-colour QCD is considered in the Baxter-Sklyanin 
representation, where the wave function is expressed as a product of Baxter functions Q{X) 
and a pseudo-vacuum state. We find n solutions of the Baxter equation for a composite state 
of n gluons with poles of rank r in the upper A semi-plane and of rank n — 1 — r in the 
lower A semi-plane (0 < r < n — 1). These solutions are related by n — 2 linear equations with 
coefficients depending on coth(7rA). The poles cancel in the wave function, bilinear combination 
of holomorphic and anti-holomorphic Baxter functions, guaranteeing its normalizability. The 
quantization of the intercepts of the corresponding Regge singularities appears as a result of 
the physical requirements that the holomorphic energies for all solutions of the Baxter equation 
are the same and the total energies, calculated around two singularities A, A* ±i, coincide. It 
results in simple properties of the zeroes of the Baxter functions. For illustration we calculate 
the parameters of the reggeon states constructed from three and four gluons. For the Odderon 
the ground state has conformal spin \m — m\ = 1 and its intercept equals unity. The ground 
state of four reggeized gluons possesses conformal spin 2 and its intercept turns out to be higher 
than that for the BFKL Pomeron. We calculate the anomalous dimensions of the corresponding 
operators for arbitrary agjuj. 



1 Introduction 

The leading logarithmic asymptotics (LLA) of scattering amplitudes in the Regge limit of high 
energies ^/s and fixed momentum transfers q = \f—i is obtained by calculating and summing 
all contributions {g^\iis)\ where g is the QCD coupling constant. In this approximation the 
BFKL Pomeron is a composite state of two reggeized gluons [|I|. The BFKL equation for the 
Pomeron wave function is closely related to the DGLAP equation for the parton distributions 
1^ . Next-to-leading corrections to its kernel were calculated in QCD ^ and in supersymmetric 
gauge theories 0. 

The asymptotic behaviour oc s-'" of scattering amplitudes is governed by the j-plane singu- 
larities of the t-channel partial waves /j (t) 
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Here the contour of integration in j is situated to the right of the leading singularity jo = 1 + 1^0 
of fj{t) (jo < 0") and the signature factor is ~ i""^ for the t-channel exchange of n-reggeized 
gluons. Its intercept ujq is proportional to the ground state energy Eq of the Schrodinger-like 
equation [jl], ||]: 

Hf = Ef ,u;o = -f^N,Eo. (2) 

The wave function / depends on the two-dimensional impact parameters - positions of the 
reggeized gluons. It is convenient to introduce the holomorphic (p^ = Xk + iyk) and anti- 
holomorphic (p^ = — iyk) coordinates and their corresponding momenta pk = and 

In multicolour QCD Nc ^ 00 the colour structure of the BFKL equation in LLA is sig- 
nificantly simplified. As a result, each reggeized gluon interacts only with its two neighbours 
i: 

1 

H = - Yl Hk,k+i ■ (3) 
^ k=i 

Note, that for three gluon composite state, describing the Odderon responsible for the high 
energy behaviour of the differences of total cross-sections for particle-particle and particle-anti- 
particle interactions, this simplification is valid for arbitrary Nc P, 0] • 

The Hamiltonian H has the properties of the Mobius invariance |^ and of the holomorphic 
separability 0: 

H=^-{h + h*), [h,h*] = 0, (4) 
where the holomorphic and anti-holomorphic Hamiltonians 

n n 

h = hk,k+i , h* = hl,k+i (5) 

k=l k=l 

are expressed in terms of the BFKL operator ^ 

hk,k+i = log(pfc) + log(pfc+i) + — log(pA;,fc+i) Pk + ^— log(pfe,fc+i) pk+i + 27. (6) 

Pk Pk+l 

Here pk,k+i = Pk — Pk+i and 7 = — is the Euler-Mascheroni constant. 

The wave function fmmi'Pii P25 ■■■■,'Pn\ Po) of the colourless composite state described by the 
operator (po) belongs to the principal series of unitary representations of the Mobius group 
For these representations the conformal weights 

m = 1/2 iz/ n/2 , m = 1/2 + iv - n/2 (7) 

are expressed in terms of the anomalous dimension 7 = 1/2 + ii/ of (po) and its integer 
conformal spin n. Furthermore, the eigenvalues of two Casimir operators and M*^ of the 
Mobius group are equal to m(m — 1) and m(m — 1), respectively. 

Owing to the holomorphic separability of the wave function has the property of holo- 
morphic factorization 0: 

/m,m(Pl> P2, •-, P^t; Po) = X! Cr,/ /m(Pl> P2, Pn] Po) f~{p*l, P*2, Pn, P*o) , (§) 



where r and / enumerate the different solutions of the Schrodinger equations in the holomorphic 
and anti-holomorphic sub-spaces: 



fm ^ f-i 



E. 



m 1 m,m ' m 



(9) 



Similarly to the case of two-dimensional conformal field theories, the coefficients Crj are obtained 
imposing single- valuedness to /mm(pi) P2; ■■■,Pn', Po) as a function of the two-dimensional vari- 
ables T^^j. 

There are two different normalization conditions for the wave function [Pl: 



r=l 



n Prvr+l / 



r=l 



r=l 



EPr f 



r=l 



(10) 



compatible with the hermiticity of H. This property is related with the fact 0, that h commutes 
with the differential operator 

A = pl2p23---PnlPlP2---Pn- (H) 

Furthermore there is a family {qr}of mutually commuting differential operators being the 
integrals of motion: 

[qr,qs]=0 , [qr,h]=0. (12) 
They can be obtained by the small-w expansion of the transfer matrix for the XXX model [|TU| 



T{u) = tr [Li{u)L2{u)...Ln{u)] = E^""'?^ 



r=0 



where the L-operators are 



Lk(u) 



u + PkoPk -Pk 

PkoPk U-pkoPk 



u 
u 



In particular is equal to A and q2 is proportional to M^. 

The transfer matrix is the trace of the monodromy matrix t[u): 

T{u) = tr [t{u)] , t{u) = Li{u)L2{u)...L„{u) . 

It can be shown that t{u) satisfies the Yang-Baxter equation: 

(n) {v) iZl {v-u) = 1:1:1 {v - u) 4 (v) 4 {u) 

where l{w) is the L-operator for the well-known Heisenberg spin chain: 

Q2M = ws:i62 + tS2S2. 



(13) 



) + (;j(-p.o i)Pk. (14) 



(15) 



(16) 



;i7) 



2 Baxter- Sklyanin representation 

Thus, the problem of finding solutions of the Schrodinger equation for the reggeized gluon 
interaction reduces to the search of a representation for the monodromy matrix satisfying the 
Yang-Baxter bilinear relations |l^ . For this purpose the algebraic Bethe Ansatz is appropriate 



11 1 . It is important, that the reggeon Hamiltonian in the multi-colour QCD coincides with the 



local Hamiltonian of the integrable Heisenberg model with the spins being the generators of the 



non-compact Mobius group SL{2,C) |jT2|, 13]. For the case of three gluons the wave function 



and intercepts of the corresponding Regge singularities were found in ref.|jTj] with the use the 
integral of motion discovered in ref.0. 

The integrals of motion and the hamiltonian for n reggeized gluons have an additional 
symmetry under the transformation W^ 



Pk Pk,k+i Pk+i , (18) 

combined with the operator transposition. This duality symmetry allows to relate the wave 
function of a composite state with the Fourier transformed wave function of (generally) another 
physical state. In particular, it gives a possibility to construct a new Odderon solution having 
the intercept exactly equal to unity . The duality symmetry can be interpreted as a symme- 



try among the states constructed from the reggeons with positive and negative signatures [p!q| . 
Indeed, the Regge trajectories for these states with gluon quantum numbers are degenerated 
in multi-colour QCD. 



In the framework of the Bethe Ansatz it is convenient to work in the conjugated space 
where the monodromy matrix is parametrized as follows. 



Here Lk{u) is given by 



t>)^z„(.)...z.(.)^(^(;) ^S)- (19) 



L,{u)=r^^f'' ]■ (20) 

V PkPko U-pkpkO I 



Now the equation for the pseudo-vacuum state 

C{u) |0)* = (21) 



has the following solution ||T3| 



|0)*=n-r- (22) 

k=l PkO 

In the total impact parameter space Jf the pseudo-vacuum wave function, 

n 1 

^'(°)(p?, pi, ...,p:;p^) = n^- (23) 

fc=i \Pm 

is an eigenfunction of the transfer matrix, 

[A{u) + Diu)] \oy = [{u-ir + iu + i r] \oy. (24) 

The pseudo-vacuum state does not belong to the principal series of unitary representations 
because it has the conformal weight m = n. 

A powerful approach to construct the physical states in the framework of the Bethe Ansatz 
is based on the use of the Baxter equation for the Baxter function Q{X) [T^, |T^. The Baxter 
equation for the n-reggeon composite states can be written as follows (see [|1^, |1^, pO|l ) 



A(") {X; fl)Q (A; m, p) = {X + 0" Q {X + t; m, p) + {X - z)" Q {X - t; m, fl) , (25) 



where A'^")(A) is the polynomial 



A(")(A; /2) = 5:(-^)VA.A"-^ fio = 2, /^i = 0, /is = m(m - 1) 

k=0 



(26) 



Here we assume in accordance with ref. that the quantities fik = <lk are real which is 
compatible with the single-valuedness condition of the wave functions [qk stand here for the 
eigenvalues of the integrals of motion]. 

The eigenfunctions of the holomorphic Schrodinger equation can be expressed through the 
Baxter function Q{X) using the Sklyanin Ansatz 



/(pi,P2, ...,p„;po) = Q(Ai; m,fI)Q(\2; m, fl)...Q(\n~i; m,/i)|0)*. (27) 
where A,, are the operator zeroes of the matrix element B{u) of the monodromy matrix: 



n-l 



Biu) = -P Ui^-Xr), P=Y.Pk. 

k=l 



r=l 



These expressions are well defined because the operators A^ and P commute |T8 



A-r , As 



Ar ) P 



0. 



(2^ 



(29) 



In ref. [0 it was assumed without any convincing arguments, that Q{X) is an entire function 
in the complex A plane. One of the purposes of our previous paper [|1^] was to find the class 
of functions to which the Baxter functions belong. For this purpose we performed an unitary 
transformation of the wave function of the composite state of n reggeized gluons from the 
coordinate representation to the Baxter-Sklyanin representation in which the operators A^ are 
diagonal [^] (see also [^). The kernel of this transformation was expressed through the 
eigenfunctions of the operators B{u) and B*{u). For the cases of the Pomeron and Odderon 
the unitary transformation was constructed in an explicit form [|l^]. As a consequence of the 



single-valuedness condition for the kernel of this transformation one obtains the quantization of 
the arguments of the Baxter functions Q{X) and Q{X*) in the holomorphic and anti-holomorphic 
sub-spaces (see [IT^ PU[): 

A^ 

X = a + i- , X* = a-i- , (30) 

where a and A^ are real and integer numbers, respectively. 

In ref.|]T^ we proposed a general method of solving the Baxter equation for the n-reggeon 
composite state. To begin with, the simplest ra-reggeon solution of this equation is searched in 
the form of a sum over the poles of orders from 1 up to n — 1 situated in the upper semi-plane 



Q("-i) (A; m,fl) = Y: 



(n-2) 



r;Tn,fi 



(A) 



t^o (A - « r) 



(31) 



where the Pr^m'^JW polynomials in A of degree n — 2. Inserting this Ansatz in the Baxter 
equation leads to recurrence relations for the polynomials P^^.^ ^^(A), which allows us to calculate 



them successively starting from -Po.mJ^(A) ||T9| 



One can normalize this solution imposing the constraint 



lim Pr 



(n-2) 



0,m,p, 



(A) 



1 



(32) 



Then, the remaining coefficients of the polynomial PqIu'^JW calculated from the condition 

lim Q("-i) (A; m, fl) ~ A-"+™ (33) 



which is a necessary condition for Q*-" (A; m, fl) to be a solution of the Baxter equation at 
A ^ oo. It is enough to require 

OO P'^"~^VA') 

lim A"-2 y ^ = . (34) 

This condition gives n — 1 linear equations allowing to calculate all coefficients of the polynomial 

<":?(A). 

The existence of the second independent solution 

Q(°) (A; m, fl) = Q^^^^ (-A; m, /?) = f , (35) 

r=0 I ^ ^'^J 



where 



= (-i)Vfc: 



is related with the invariance of the Baxter equation under the simultaneous transformations 

A — ^ —A , /i — > /i^ . 

One can verify |]19| that 

Q* (-A; m, /T'*) = Q {X*; m, /i 

It turns out |Il9|, that there is a set of the Baxter functions Q*^*^ {t = 0, 1, ...,n — 1) having 
poles simultaneously in the upper and lower half-A planes. 



(A; m,/2) = y 



r=0 



{\-iry (-A-ir)"-i- 



where the polynomials P^*'^^^ and i^^^"^^^*) are fixed by the recurrence relations following from 
the Baxter equation and from the condition that the above Baxter functions decrease at infinity 
more rapidly than A"""*"^. These solutions Q*-*^ are linear combinations of Q'^'^~^^ (A; m, /T) and 
A; m, with the coefficients depending on coth(7rA) ||19| . 
Using all these functions in the holomorphic and anti-holomorphic spaces one can construct 
the normalizable Baxter function A ^ in the space (cr, A^) without poles at cr = 

Q^,^r^{J) =Y.Ct,iQ'^'^ (A; m,/i) (a*; m, /?) (36) 

t,i 

by adjusting for this purpose the coefficients Ct^i. 

Another problem in the Baxter approach is the calculation of the energy E^~, because the 
expression suggested in ref. |jl3| leads to an infinite result for meromorphic Baxter functions. 
This problem was solved in ref.[|T^ by an unitary transformation of the Hamiltonian to the 
Baxter- Sklyanin representation. 



In the region, where the gluon momenta are strongly ordered in their values, 



\Pn\ « bn-l| << ••• << bl| = 1 

the unitary transformation of the wave function \E'^ ~ to the Baxter-Sklyanin representation is 
significantly simplified [jl9| 



k=i y-°° Nk=-oo 



On the other hand, in this region 

2 1 

^m,mbl.--- >Pn) ~C„ |p„| In- ^, (38) 

\Pn\ 

where Cn is a constant [|l^. Therefore \E'^ ~(Ai, . . . , A „_i) has the first order poles at A„,_i = i 
and Xn^i = i for Nn-i = and the pole singularities at A^ = A^ = (fc = 1,2, ...,n — 
2). The action of the Hamiltonian H on the function ^E'^- near these singularities after the 
unitary transformation to the A-space is drastically simplified [|19|. Moreover, using the Sklyanin 
factorized expression for the wave function ~ (Ai, . . . , A n-i) we express the energy in terms 



of the behaviour of Q(A„_i) near the pole at A„_i = z, A*_]^ = i and obtain [|T^] 

d d 



E = i lim -— — — In 

A,A-^i d\ dX* 



{x-^r-\x*-^r-'\xrQ^,^,p{'X)] . (39) 



It is important, that the consideration of the kinematical region with an opposite ordering of 
the gluon momenta pk leads to an analogous expression for E in terms of the behaviour of 
Qm,rn,pi ( ) ^^ar A, A* = —i. The coincidence of these energies imposes a constraint on the 
spectrum of the integrals of motion {k = 3,4, ...,n). 

Since the function Q^n - p; A ^ is a bilinear combination of the independent Baxter functions 
Q^^\X) and Q^''\X*) {t,l = l,2,...,n), the holomorphic (anti-holomorphic) energies for all 



solutions should be the same 



i lim In 

A-i dX 



This leads to the quantization of the integrals of motion • 

Note, that this condition can be obtained also as one of the consequences of the absence 
of poles in Qm,m,fi ( ) at 0" = for |A^| > 0. Because the products of poles (A — ir)~^ and 
(—A* — ir'Y' appear in the corresponding bilinear combinations, to cancel all singularities one 
should expand the functions Q^^\ Q^^^ in Laurent series near these poles and keep all relevant 
terms. The coefficients of the Laurent series satisfy recurrence relations similar to those for the 
residues of the poles. However, it turns out that they are not proportional to the residues of 
the poles.. Even a special choice of the values for the integrals of motion qr does not lead to 
such proportionality because these recurrence relations are different at = 2 ||19|| . 



We discuss below the solution of the Baxter equation by the method proposed in ref. 
Note, that for n = 3 this method gives results in full agreement with those obtained by other 



approaches flj, |15[. We find in the present paper the wave functions and intercepts for the 



quarteton (four reggeons state). 



3 Meromorphic solutions of the Baxter equation 

Let us rewrite the Baxter equation in a real form introducing the new variable x = —iX, 

n{x, pi) Q{x, jl) = {x + 1)" Q{x + 1, /i) + (x - 1)" g(x - 1, /i) , (41) 

where 



Q{x,fl) = J2 (-l)Vfc^' 

k=0 



n—k 



(42) 



and 



/io = 2 , /ii = , /i2 = m{m - 1) , 

assuming that the eigenvalues of the integrals of motion /i^ {k > 2) are real numbers. 

Then, the solution of the Baxter equation possessing poles only at a; = / = 0, 1, 2, ... can be 
written as follows, 



Q 



(n-l). 



X. 



=0 I 



+ 



(x - /)"-! {x - ly 



+ ... + 



X — I 



(43) 



where the residues a/ (/I), bi{ll), zi{jl) satisfy recurrence relations which can be easily obtained 
from the Baxter equation and its derivatives in the limit x I. Using these relations we can 
express all these residues in terms of ao{fI), 6o(/^), ZQ^jl). In addition, imposing the Baxter 
equation at x — *• cxd. 



lim x" g("-^^(a;,/2) = , s = 1, 2, . . . , n - 2 ^ 



(44) 



fixes the parameters bo{fI),..., ZQ{fI) in terms of aQ{fl). We then require the normalization 
condition 

ao(/2) = 1 , 

without losing generality. In this normalization the holomorphic energy is (see |T9| ) 

bl , , fin-l 

e„ = \-n = bo . 

Thus, the solution Q^^~^\x,jl) is uniquely defined and can be explicitly constructed. 
Let us further introduce a set of the auxiliary functions for r = l,2, ...,n — 1 



1=0 



ai{fl) hi{p) 



{x-iy {x~iy-^ 



X — I 



(45) 



where the coefficients ai,...,gi satisfy the same recurrent relations as ai,...,zi but with other 
initial conditions 

ao = l,&o = -=^o = 0. (46) 

Note, that all functions fr{x, p) are expressed in terms of a subset of pole residues ai, ...,zi for 
fn-i{x,p). 

Now we write the Baxter function Q^^~^\x, p) as a linear combination of fr{x, p) 



Q^''-'\x,P) = Y.Cr{P)fr{x,P). 



n-l 



(47) 



r=l 



In order to impose the asymptotic condition (^) on Q^" ^\x, fl), we use the binomial series 
for the terms in eg. (|^) , 

1 _ ^ /" (j + ri- 1)! _ ^ /'^-J' 

- ly ~ h> ^! (j - 1)' Mj - 1)! ■ 

We find a set of n — 2 hnear equations on the coefficients Cr{p)'- 

n—l 

Y,GkAf^)Crifl) = , A; = l,2,...,n-2 = 1 , 

r=l 

where 



1=0 



(4J 



Obviously, the symmetric solution 

g(o)(a;,/i) =g("-i)(-x,/?) , where /i^. = , (49) 

can be constructed in a similar way. It has poles at x = — / (/ = 0, 1, ...). 

But there are other 'minimal' solutions Q^^\x, p) (t = 1, 2, n — 2) of the Baxter equation 
having t-order poles at positive integer x and (n — l — i(:)-order poles at negative integer x [|19| 

t n-l-t 

QW(x,/2) = 5:cW(/2)/,(x,/i)+/3W(/i) Y: Cl-~'-'\fi) U-xJ^) , (50) 

r=l r=l 

where the meromorphic functions fr{x,fl) were defined above. Such form of the solution is 
related by the invariance of the Baxter equation under the substitution 

The coefficients CW(/i) , C^"-i-*)(/r^) and (3\pj are obtained imposing the asymptotic va- 
lidity of the Baxter equation, 

Yim x^ Q^\x,Il) = Q , A; = l,2,...,n-2 . 

X — >oo 

This leads to a system of n — 2 linear equations 

t n-l-t 
r=l r=l 

A; = l,2,...,n-2 , (51) 

where the matrix elements Gk,r{j^) are defined by eq.(^) and we normalize Q^^\x, pi) by choos- 
ing 

Cf\fl) = Ct-i-t\fl) = l. (52) 
Moreover, one can obtain the following relations using the symmetry of the Baxter equation, 
gW(x,/2) = /3M(/i) g("-i-')(-x,/?) , /3M(/I) /3("-^--)(/?) = 1 , /?(°)(/I) = 1 . (53) 



It is important to notice that three subsequent solutions Q*^^^ for r = 1,2, ...,n — 2 are 
hnearly related by 



S^'^^fl) + 71 cot(7rx)l Q^''\x,fl) = g(''+^)(x,/i) + a^''\fl) Q^''-^\x,Jl) . 



(54) 



Indeed, the left and right-hand sides satisfy the Baxter equation everywhere including x 
oo and have the same singularities. Due to the uniqueness of the 'minimal' solutions the 
quantity cot(7rx) Q^'^\x, fx) can be expressed as a linear combination of Q^'^~^\x, pi), Q^'^\x, fx) 
and Q^'^'^^^x , fx) . Furthermore, the coefficient in front of Q^^~^^\x, fl) is chosen to be 1 taking 
into account our normalization of Q^^\x,fx). In this normalization we obtain 



/?M(/2) = (-irn«^*H/^)- 



i=l 



The following relations hold for the coefficients a^^\fi) and 5^^\fl) 

S^'\fl) = -S^-~^~-'\i?) , a^'\fl) a^''-^-'\i?) = 1 . 



(55) 



(56) 



According to the above recurrence relations the functions Q^^\x, fx) (r = 1 , 2, . . . , n — 2) can 
be expressed as linear combinations of Q^^~^\x,fx) and Q^^\x,fx) with the coefficients being 
periodic functions of x |T^: 



D 



S{fl), d{fl) , IT COt{TTx) Q^^\x,fx) 



(57) 



(r) 



6{fl), d{fl), IT cot{TTx) Q^^\x, fl) + dI^Ii 6 {fl), a{fl), TT COt{TTx) Q''" ^\x,fl) 



compatible with the relation 



The factor D 



6, a, IT cotfvrx) 



can be written as the determinant of the matrix A 



D 



6, a, TT cot(7ra:;) 



A 



6, a, vr cot(7rx) 



(5^ 



The matrix A takes the form, 



A 



6, a, n cot(7rx) = vr cot(7rx) / + A 



6, a 



and 



A 



6, a 



( (5« -1 

^(3) 



V 







/ 1 
1 

V 



\ 





\ 


1 ) 



(59) 



(60) 



... y 



In an analogous way d'q^ 6, a, 7rcot(7rx) and -D^^i^ 6, a, n cot{'Kx) can be expressed in terms 

of the determinants of the matrices Ag*^^ and A-^^-i of (lower) rank n — 3 obtained from A by 
removing the column r and the first or the last line, respectively: 



(r) 



S,a,7i cot{TTx) = (— l)*" a*-^-* 



A 



(r) 



D^^ll 6,d,TT COt{TTx) = ( — 1) 



A, 



(r) 
n-1 



-1 r\ 



and 



A, 



(r) 



/ ^(2) +v]-cot(7rx) 





5(^+1) +7rcot( 



TTX] 









V 







/ + TTCotfvrx) -1 

+ TT cot(7rx) 



V 

Let us write these combinations in the form 



^('^-s) + 7rcot(7ra;) j 



\ 


-1 ... 



g('')(x, jl) = a(")[/i,7rcot(7rx)] g(")[x, /i] + vr cot(7rx)] g(°) [x, jl] 

for n + 1 reggeons and 

Q^'\x, Jl) = a'(")[/i,7rcot(7rx)] g^^+^H^^, fi) + b'^''^[fl,7T cot{7Tx)] Q^^\x, fl) 
for n + 2 reggeons. Then, the following relations are valid 

a'*-^^ [/?, 7r cot(7ra;)] = a'-'"-' [/I, tt cot(7rx)] a"-"-* [/I, tt cot(7rx)] , 

[/2,7r cot (ttx)] = 6('')[/I,7rcot(7rx)] + a^''\fl, n cot {nx)] 6'^"^ [/I, tt cot(7rx)] 



where 



a'^"U/i, vrcotfvrx)] = ^, . ; — , . . — -r-— , 

^ ^ 5(")(/2) + 7rcot(7rx) - a(")(/i)a("-i)(/i) 

6'(")[/i,7rcot(7rx)] = 6^"-^) [/T, vr cot(7rx)] a'^") [/?, tt cot (vrx)] 



(61) 



(62) 



Notice that the linear relations (|5^) among g*^''^(x, /I), g 



X. 



/2) andg('-i)(x, fl) are similar 
to the recurrence relations for the orthogonal polynomials Pr{z) if we substitute tt cot(7ra;) by 
a variable z. The Baxter functions also belong to an orthonormalized set of functions. 

We use below the formulae of this section for the numerical calculations of the intercepts of 
the composite states constructed from reggeized gluons. 



4 Spectrum of eigenvalues of integrals of motion 

To quantize the integrals of motion one should impose all physical constraints on the corre- 
sponding eigenfunctions. In the impact parameter representation p the wave function should 
be normalized in a conformally - invariant way [|, ^ and the Hamiltonian H in the space of 
these functions should be hermitian. The integrals of motion {ql) are differential operators in 
the holomorphic (antiholomorphic) subspaces [jlO[- The wave function in the p-representation 



is constructed in a bilinear form from all eigenfunctions of the integrals of motion similar to the 
case of two-dimensional conformal field theories in such a way that its single-valuedness holds. 



The single-valuedness condition is an important constraint on the coefficients of the bihnear 
form and on the possible values of the integrals of motion. But there are other constraints 
which are significantly simpler. 

In the case of the Pomeron one does not have extra integrals of motion apart from two 
Casimir operators with their eigenvalues depending on the conformal weights m and m. For 
the principal series of the unitary representations of the Mobius group the weights depend on 
the real parameter v and the integer conformal spin n 0. 

For three particles there is an additional integral of motion with its eigenvalue parametrized 
as gs = The problem is to find the region to which the parameter belongs. For hermitian 
operators the eigenvalues are real, but the integrals of motion in holomorphic and anti- 

holomorphic subspaces are not hermitian. With the use of conformal invariance the integral of 
motion ga for the Odderon can be written as the third order differential operator f^, 15 



z{\- z) z{\ - z)p^ + %{2 - m)(l - 2z)'p^ + (2 - m)(l - m)p 



acting on the functions depending on the anharmonic ratio 

P12P30 
z = . 

P10P32 



Here P = i-§^ and is the duality operator |[T^ 



am = z{l - z)p 



l+m 



(63) 



(64) 



(65) 



transforming a state with the conformal weight m into a state with conformal weight 1 — m 



Ir 



>l-m 



k- 



-mfm ■ 



Here Im and li-m are the eigenvalues of the corresponding transformations depending on the 
norms of the functions 0m- The product of these numbers does not depend on the normalization 
and is related with the eigenvalue by 



93 = ^P • 



The Odderon wave function 
holomorphic functions 



z) is the bilinear combination of holomorphic and anti- 



ik 



m,q[i \ 



'm,-ij3 \ 



(66) 



where the symmetrization under the transformation ^3 —q^ is needed for the Bose symmetry 
of the wave function. The coefficients Cik are calculated from the single-valuedness of 4>mm,i^- 
The sum is performed over three independent eigenfunctions of the integrals of motion q^, q^ in 



the holomorphic and anti-holomorphic spaces. It was argued in ref. that the eigenvalues q^ 
are purely imaginary as a consequence of single-valuedness. We give below additional arguments 
supporting this conclusion [at least for real values of m(m — 1)]. 

To begin with, let us perform the simultaneous interchange of {z, m) and {z*,m) in 4>mmi^)- 
The wave function should remain the same if there is no accidental degeneracy. But this can 
be true only if ^3 satisfies one of the two relations ^5 



93 = ±93 > 

-1 n 



(67) 



which is in agreement with the Janik-Wosiek result qs = —q^ [|T^|. An analogous argument for 
the general n-reggeons case will lead to the guess that fl* = fl ot fl* = fl'^ . 

Furthermore, although the holomorphic hamiltonian h is not a single-valued operator (see 
for example the discussion of the separability properties of the BFKL hamiltonian in ref. ||19||), 
its ambiguity seems to be simple and related with the possibihty to add to it some periodic 
functions of d/d{\npk) cancelling in the total hamiltonian H. For example, in the case of the 
Odderon one can write the holomorphic hamiltonian in the normal form 



h = - \n{z) + ^(1 - P) + ^(-P) + ^(m - P) - 3^(1) + M^) 



k=l 



where z is the anharmonic ratio, P = d/d{ln z) and the explicit formulae for fk{P) are given in 
ref. [15], where h was defined with an extra factor 2 in comparison with the present paper and 
with ref. [11^. The above expression for h is equivalent to the expressions obtained by adding 



to it one of two terms 



-271 cot(7rP) , vr [cot(7r(m - P)) - cot(7rP)] 



cancelling poles of the \l/ -functions. Apart from such ambiguities, h can be considered for 
real m(m — 1) as a symmetric operator h = acting on functions z*) depending on real 
variables z, z* with the norm 



dz* 



oo Z*{1 — Z*) 



dz 



zl 



The integration over the real values of z and z* can be obtained after the anti-Wick rotation 
of the vector z to the Minkowski space. Therefore the eigenvalues e of h should be real (for 
the case when m{m — 1) is real). On the other hand, e is a real function of the eigenvalue /i^ 
(apart from constant imaginary contributions cancelling in the total energy E). For example, 
we have for large /i the expansion of e in 



15 



e = ln(/i) + 37 



A + _ 1/2)2 _ _ ^2)4 

448 120^ ' ^ 12^ ' \ 




Therefore /i is real or pure imaginary for real m(m — 1) (at least for large |yu|). 

Because the operator q^ is ant i- symmetric under the permutation of the coordinates pk and 
pi, its matrix element between the symmetric wave functions is zero, but the matrix element of 
g| can be written as follows 



+00 



dz* 



00 Z*{1 — z*) 



dz 



zl 



gangs') 



and therefore the eigenvalues p of —iqs are again real or pure imaginary for real m(m — 1). 

Let us consider now the expression for the total energy of the composite state of n-reggeons 
in terms of the holomorphic and anti-holomorphic energies 



,{p) + e~{r) 



valid for the wave function 



satisfying the Schrodinger equation in the Baxter-Sklyanin 



representation in the limit X,X* i [|T^. We can obtain the analogous expression 



P. 



emip") + e~{p*) 



-1 o 



by taking instead another limit A, A* —i. These two expressions for energies were derived 
from the Schrodinger equation with the hermitian hamiltonian [19|. Therefore they should 
coincide for its eigenfunction 0„ ~. This is possible only if the following property is fulfilled for 
the quantized values of fl: 

em(/i) + e~(^'*) = e„(/i') + e~(/I*) , 
The solution of this equation for real m(m — 1) is 



/i = /i, 



(68) 



19], that the eigenvalues /i 



or jl* = jl^, if again there is no accidental degeneracy. 

The arguments given in this section support our assumption 
are real for real m(m — 1). For non- vanishing conformal spins n the values of e and /I can be 
obtained by analytic continuation of the anomalous dimensions of the corresponding high-twist 
operators with n = to continuous values of the Lorentz spin j (see sections 6-8 below). 

In our opinion the complex values for /i found in do not correspond to eigenfunctions 
of the Schrodinger equation for the Odderon and they were obtained due to an incorrect quan- 
tization procedure. 



5 Zeroes of the Baxter function and the quantization of 
the integrals of motion 

The zeroes of the Baxter function are very important, because their position Afc is fixed by 
the Bethe equations and their knowledge gives a possibility to write the wave function of the 
composite states in the framework of the Bethe Ansatz as an (infinite) product of the differential 
operators -B(Afc) applied to the pseudo-vacuum state. The positions of the zeroes of Q{X) for 
the Pomeron wave function was investigated in our previous paper For the case of the 



composite states constructed from n > 2 reggeons the number of the 'minimal' solutions Q*-*^^ of 
the Baxter equation is n and their linear combinations have rather complicated sets of zeroes. 

However, certain linear combinations of functions Q^'^\x, fx) have zeroes which are situated 
at equidistant points Xk = xo + k, k E Z. Indeed, according to relations ( ^i)) among the Baxter 
functions Q^'^\x^ jl), Q^'^~^^\x, jl) and Q^'^~^\x, jl) such situation happens for solutions of the 
equation 

g("+i) (x, /i) + (/i) g('--i) (x, /i) = , (69) 



where r = 1,2, ... ,n — 2. Some of the solutions of eq.(|5UD coincide with the zeroes of the 
function Q^'^\x, jl) while other roots are situated at the equidistant points 



(/I) = /c arccot 



TT 



5M(/i) 



TT 



keZ. (70) 



For the eigenstates of the Hamiltonian the parameters fl are quantized in accordance with 



the physical requirement |19| that the holomorphic energies for different Baxter functions are 



the same. It is equivalent to the condition, that all parameters 5^'^\jl) vanish 

5W(/2) = 0, (71) 



-1 A 



because otherwise, the energies for the solutions Q^^\x,fl) and Q''^^^\x, fi) would not coincide. 
Eq.([fOD then implies that for quantized /i the above linear combination of Q^'^^^\x,fl) and 
Q{r-^)(^x, pi) has a sequence of zeroes at the points 

Xfc = /c + ^ , ke Z. 

Let us consider now the superpositions of the Baxter functions Q^^~^\x, p) and Q^^\x,p) 
with their poles situated only at positive and negative integer points, respectively. It is obvious 
from the previous section that there are n — 2 different linear combinations of these functions 



<l>W(x,/2) = Q^''-'\x,P) + c^'\p)Q^'\x,P) 
having equidistant zeroes at the points x = Xk where (see (58)) 



D 



S{p), d{p) , 71 COt^Tlx] 



The last equation has n — 2 different solutions 

Z^^\P) = TT cot TTX^^^p) 



t= 1,2, 



.n 



(72) 



(73) 



(74) 



and for each solution there is a linear combination of the Baxter functions Q^" ^\x,p) and 
Q^^\x,p) with the relative coefficient 



cW(/2) 



(r) 



5{p),a{p) ,z^^\p) 



D 



(r) 
n-1 



6{p),d{p),z(^\p) 



(75) 



Note, that c^^\p) does not depend on the parameter r, if there is no accidental degeneracy. 
In the case, when 

5^''\p)=0 (76) 

for all r = 1, 2, n — 2, which corresponds to the quantization of the integrals of motion p, the 
determinant D S{p), d{p), 7rcot(7rx) is an even (odd) function of its argument 7rcot(7ra;) for 

even (odd) n. Therefore, the equidistant zeroes of two different functions ^^^\x,p) (or zeroes 
of the same function ^^^\x,p)) have opposite signs (modulo an integer number): 



x^'\p) 



-x^^-'-'\P) 



(77) 



For odd n one function has a sequence of zeroes at 



In a general case of n reggeized gluons we can calculate the important coefficients a^'^\p) 
in the recurrent relations and the quantized values of p only by finding Q^"'\ Q^^^ and all their 
linear combinations with constant coefficients which have the equidistant zeroes with above 
properties. Let us consider several examples: 

For n = 3 we have only one function with equidistant zeroes Xk = k + 1/2 



Q{x) = Q(2) (a;, ^) + (^) g(o) (^.^ ^) 



For n = A there are two functions, 

with zeroes at tt cot ttx = ±^/a2, respectively. For n = 5 there are also two functions with 
equidistant zeroes: 

Qi = Q(^) + a(^) g(°) , = Q^'^ - a^^') Q(°) 

with the equidistant zeroes at vrcotTra; = ±V a^^) + and vrcotvrx = 0, respectively. In all 
these cases, one can calculate all a^"^^ and the quantized values of /2 by finding the corresponding 
combinations of Q*-"^ and Q^^^ with the equidistant zeroes. 

6 Anomalous dimensions of quasi-partonic operators 

The Q^-dependence of the inclusive probabilities ni{x, In Q"^) to have a parton i with momentum 
fraction x inside a hadron with large momentum |"p^| — > oo can be found from the DGLAP 
evolution equation 0. The eigenvalues of its integral kernels describing the inclusive parton 
transitions i ^ k coincide with the matrix elements 7j^*(a) of the anomalous dimension matrix 
for twist-2 operators with the Lorentz spins j = 2, 3, ... . These operators are bilinear in the 
gluon {i = g) or quark {i = q) fields. For example, the twist-2 gluon operator with the Lorentz 
spin j can be written as, 

01 = n^'^n^\..n^^ tr G,^,D^,D^.^...D^^_,Gp^. , (78) 

where = + gV^, and Gp^ = ^ [Dp,Dfj] are the covariant derivative, the gluon field 
and the field tensor, respectively. 

The symmetric traceless tensor O^^^j...^^ is multiplied by the light-cone vectors n^^ 

I 2 n 2 n 2 ,o2 Q 

n^, = Qp + xpp, n^ = 0, ~ , qp = -Q , ^ = ^ ' 

where and in the deep-inelastic ep scattering are the momenta of the initial proton and 
virtual photon, respectively. 

The matrix elements of the operators between the hadron states are renormalized as 
functions of the growing ultraviolet cut-off Q^. For example, in the case of the pure Yang- Mills 
theory with the gauge group SU{Nc) we have 



p 



\ ( rQ^ \ Att 11 2 

~ exp / j,{a,{Q"))d\nQ" , a^iQ') ^ — , P2 = -^N, - i^rif 
' VQl J /?2ln^ 3 3 



where — 200 Mev is the QCD parameter and the anomalous dimension 7j(a) can be calcu- 
lated perturbatively as 

7» = £cf'(^)\ (79) 

k=l ^ 



TT 



In particular, to the lowest order 

J j + 1 j + 2 12 



-1 r> 



The gluon anomalous dimension is singular in the non-physical point u = j — 1 ^ 0. In this 
limit one can calculate it to all orders of perturbation theory |Q 

7.. = ^-^" 1 ^ +••• 80 

from the eigenvalue of the kernel of the BFKL equation in LLA at n = 0: 

a N 

UJBFKL = ^[2^(l)-^(7)-*(l-7)] • (81) 

TT 

Notice that the coefficients of and exactly vanish in eq.(^). Here \E'"(1) = 

-2C(3) = -2.4041138. . .. 

Indeed, the Green function satisfying the inhomogeneous BFKL equation can be written as 
follows [§] 

< 't>{n)m)m)4>(0) >=i: c(.,n) /^^. ^^"■''"^■'"^'f-'f"'"^'' . (82) 

where uj^{n^ v) is the eigenvalue of the BFKL kernel and C{v,n) is fixed by the completeness 
condition for its eigenfunctions 

Ey^n{plQ,p20) =< O|0(pi)0(p2)C'^,n(Po)|O > = 

\pioP2oJ \ptoP2oJ 2 2' 2 2 ^ ^ 

In the limit \pi'2'\ one can perform a Wilson expansion for the product of the fields 
(f){pi') and 4>{p2')- In this case the integral over can be calculated by extracting the factor 
Eu,n{pio,P2o) from the integral at the point po = Pi' and by using the relation 

Pi'2'V I |2r -p m + m 

Moreover at |pi'2'| one can then shift the integration contour into the lower half of the 
i/-plane up to the first pole of [u — uj^{n, z/)]~^ having Imu < P, |23| : 



lim < 0(p?) 0(p|) 0(pr) 0(p?) >^ 
|Pi'2'Ho 



r 

where e ~ ag/uj and 



e+ioo 



C(7,n) d'j 
uo — uo^in^ v) 



P12 Pl'2' 



Pll'P22' 



2(1-7) 



7 = ^ - zz/ , (84) 



I PlV P 22' Pl'2 PV 2' 



V Pl2Pl'2'Pll'P22' 

Note, that in accordance with the fact that the Green function includes the external gluon 
propagators, the scattering amplitude behaves in the momentum space as for large 

gluon virtualities \k\^ . 



TTLJ 



In particular, for n = and 7-^0 one can obtain the above expansion ( pUD in powers of 

for the position 7^ of the pole (7 — 7aj)~^- It is important that for real uj the pole is situated 

on the real axis. Therefore, the description of the corresponding state in the framework of the 

Mobius group approach requires the exceptional series of unitary representations p5| (contrary 

to the Regge kinematics, where the principal series is used). For the exceptional series the 

anomalous dimension _ 

m + m 
7 = 1 

is real. 

One can calculate from the BFKL equation also anomalous dimensions of higher twist 
operators by solving the eigenvalue equation near other singular points 7 = —k {k = 1,2, ...) 
or by including in it a dependence from the conformal spin |n|, which also leads to a shift of 
the pole position to 7 = —\n\/2 (see 

But a more important problem is the calculation of the anomalous dimensions for the so- 



called quasi-partonic operators |2|] constructed from several gluonic fields. Indeed, the contri- 
bution of these operators at j ^ 1 is responsible for the unitarization of structure functions at 
high energies. The simplest operator of such type is the product of the twist-2 gluon operators 



p 



0' = E0':.., J = J2jr=P + UJ, U = Y.UJr. (85) 



r=l r=l r=l 



In the limit A^^^ ^ 00 this operator is multiplicatively renormalized ||26[ and its dimension 



is the sum of dimensions of its factors (including their anomalous dimensions 7(1-1;^)) 

p 

r = p -7, 7 = Y^^i^r) , (86) 

r=l 

where in the expression for the total dimension F we neglected the small contribution ^. 

To investigate the multi-Pomeron configuration, let us write the hamiltonian describing the 
corresponding composite state in LLA as a sum of pairwise BFKL hamiltonians 

p 

r=l 

neglecting the relative Pomeron interactions. In this case its eigenfunction as a product of the 
Pomeron wave functions: 



^(Pai,Pfei,Poi; •••;Paj,,Pbj,,PoJ = n 



P / \ / * * \ 

ParOrPbrOr \ / ParOrPbr-Or \ 



r=l 



Parbr / \ Parbr 



where pa^, ph^ and po,. are the coordinates of gluons a^, K and Pomerons Or, respectively, and 
the quantities mr and are their conformal weights. The wave function ^ belongs to a 
reducible representation of the Mobius group and can be expanded in a sum of irreducible 



representations with the use of the Clebsch- Gordon coefficients C"*'~ ~ [OS 

^ m-s_mi;...\mrmr ^ — ' 



\{o^^~Xyor) = E /^'poC™;^^. .^^-^(p'op p'oa— A; /o) <^'™'"'(Po)- 



The contribution to the scattering amphtude in the coordinate space from each irreducible 
component can be written as follows 



n 



r=l 



oo 



Parbr 



.ParOrPbrOr 



^PbrOrParOr 



mimi;...;m,rmr 



POp! Po) 



In the Regge regime rrir and rrir belong to the principal series of the unitary representations 



of the Mobious group and therefore m and m also belong to the same series ||2^. It means, 
that the position ujq of the t-channel partial wave singularity related to the asymptotics of the 
cross-section at ~ equals c^o = V^bfkl- 

In the deep-inelastic regime the essential intervals are small |poOrl ~ 1/Q ^ IPa^bJ- From 
dimensional considerations we obtain at large \Q\ after integration over the essential region of 
Po^ a power asymptotics of the scattering amplitude 



r=l 



^ / dUr 



Par br 



Parbr 



.ParOPbrO, 



.PbrOParO, 



mimi;...;mpmp 



\Q\ 



-2(p-7) ^inip 



i7) 



where e"^ = JQ/Q*, n = | I]r=i("^r ~ ^r) and 



7 = p - - ^(m^ + m^) (88) 

is the anomalous dimension of a composite operator which turns out to be real and small for 
small g'^/uj. This expression for 7 is in agreement with the known result ||2^, that in the 
product of the exceptional representations with real 7.^ = 1 — > there is a continuous 

spectrum of the unitary representations of the principal series and only one representation from 
the exceptional series having 



7 



r=l 



^9) 



Since u^n,^) has a negative second derivative uj'!^^{n,v) at < 7 < 1, we obtain after the 
integration over v^. with the use of the saddle-point method the following result for A in the 
Bjorken regime 

p 



Parbr 



Par br 



.ParOPbrOj \PbrOParO, 

.(-^)(c./p). 



\Q\ 



-2(p-7) 



where 

■y = p uj 

and uj^^^\uj) is the inverse function to u; = lobfkl{i)- 

In particular, for very large corresponding to 7 — * and n = we have 



(90) 
(91) 
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dsNc 2 
V 

TXUJ 



(92) 



It is possible to calculate also a correction of the relative order A^~^ to this expression . 

In accordance with the topics of this paper let us consider now the high energy asymptotics 
of the irreducible Feynman diagrams in which each of n reggeized gluons at A^^^ — > 00 interacts 



only with two neighbours. In the Born approximation the corresponding Green function is the 
product of 
follows (cf. 



product of free gluon Green functions Y[r=i \Pr ~ Prf- In LLA — ~ 1 it can be written as 



< (j){-pt)...(f){p^)(f){p[)...(f){p'J >= 



'm,m;/J3,., 



(93) 

where u{m, m; /is, //„) ~ as are the eigenvalues of H depending on the quantized integrals of 
motion /is,...,/^^ and ^m,m-^^3,...,i,„i'P'i^ p"™; p'o) are the corresponding eigenfunctions. 

From the Mobius invariance of the Schrodinger equation we obtain 

m m 

VI/ ~ ^ . ^ \ - . ( Pl2P23-Pnl \ " / P*12P*23-Pnl \ " . ^ ^ 

^m,m;/X3,/i2,.--,MnV 1' ri' oJ ~ I 2 „2 „2 ) I „*2 „*2 ) /m,ml^l5---5 ^ n-2) , 

\Pl0P20---Pn0 / \Pl0P20---Pn0 / 

(94) 

where af ,. (x,. and x*) are independent anharmonic ratios of pk and pi for A; = 0, 1, ...,n. 
In the Bjorken region, where 

\p'r- Ps\ ~ Q'^ < |Pr -Ps| , 

the essential integration domain is \po — Pr\ ~ 1/Q and therefore from dimensional considera- 
tions 



< 0(pt)...0(2) E r ^^EE-E g».,f?.;.3,... ^-'--^-'-.^^^i-' . (95) 

- 71;^ Tit u;-^(m,m;/i3,...,/i„) 



It means, that in this limit the contour of integration over u should be shifted to a lower half 
of the complex plane up to the first pole of [uj — uj{m, m; fi^, /i2, /x„)]~\ 

For small coupling constants as this singularity for n reggeized gluons is situated near the 
pole of u{m,m] /is, ...,//„). The position of the leading pole is 



2 2 ' to \l CO 



asN, 



(96) 



This relation is in agreement with the above estimate of the anomalous dimension for the 
diagrams with the t-channel exchange of p BFKL Pomerons because here n = 2p and c*^^^-* = p"^. 
Moreover, it can be obtained from the solution of the equation for matrix elements of quasi- 
partonic operators written with a double- logarithmic accuracy in ref. Let us derive a 

similar equation starting from the Schrodinger equation for the composite states of reggeized 
gluons in the multi-colour QCD. 

In the case of two reggeized gluons in the momentum space we obtain for \p[\ ~ \p2\ ^ 
\pi \ , \p2\ the equation 

^^(Pl, P2) = 1 ,2, ,2 / , ,2, ,,d\p,\ / — -^(pi,-pi). 



Note, that the Bethe-Salpeter amplitude \I'("??i, "p^2) contains the gluon propagators \pi\ ^ , \p: 
The solution of this equation is 



2 



r\r\ 



In the case of n reggeized gluons in the multi-colour QCD apart from the product of the 
propagators l\r=i \Pr\~'^ due to the gauge invariance one can extract from the amplitude also 
the momenta p^'' for each gluon 

^( P^i, 'P'2, P'n) ~ n ri2 /m,M2,...,Mn • (97) 

r=l \Pr\ 

The factor nr=iPr' / \Prf leads after the Fourier transformation to a singularity in near 
n/2. The function /^i,^2v,Mn ^ tensor depending on = In p^. The Schrodinger equation 
for this function takes the form 



The above equation gives a possibility to calculate the anomalous dimension 7 = 7(0;) = 
Cnas/uj. A similar equation is discussed in ref. [^. In particular, for the Odderon (n = 3) 
it turns out that C3 = according to an unpublished result of M. Ryskin and A. Shuvaev. 
We confirm this result below [see eq. (|117|) ] by solving the Baxter equation and finding a pole 
singularity near of ^2±ni = 2 (instead of 3/2 as it could be expected from above considerations). 
For n = 4 in accordance with the general formula "'"^"^ ~ f ~ 7^"'' a pole singularity 



near Hi±™ = 2 as shown in eq. (|130|) . Moreover, similar to the case of the BFKL Pomeron 



the anomalous dimensions 73 and 74 are calculated for arbitrary a/ 00, which is important for 
finding multi-reggeon contributions to the deep-inelastic processes at small Bjorken's variable 
X. We plot 7ruj/[a Nd as a function of m = m for the Odderon in fig.|l| and for the four-reggeon 
state in fig. ^ It is important, that for m = the curves allow to calculate the intercepts of 
some states with conformal spin k. 



7 Solutions of the Baxter equation for the Odderon 



We construct explicit solutions of the Baxter equation for the Odderon [n = 3) following the 
general method presented in sec. 3. 

The Baxter equation for the Odderon takes the real form ( ^T] ) 



where 



(x; m, fi) = + m(m — l)a; + fi Q {x; m, fi) 
-{x + 1)^ Q{x + 1; m, fi) - {x - 1)^ Q {x - 1; m, ^) = 

qs = ifi , Im (fi) = 



(99) 
(100) 



and jj is assumed to be real, which is compatible with the single-valuedness condition for the 
Odderon wave function in coordinate space IT^ . 



The auxiliary functions (^) for the Odderon take the form 

oo 

/2 (x; m, fi) = J2 

1=0 

oo 

/i (x; m, /i) = ^ 



ai{m, li) ^ bi{m, /i) 



{x-iy 



X — I 



1=0 



X — I 



(101) 



Imposing the Baxter equation i?3 = at the poles x = r yields the following recurrence relations 
(cf. [|T9l) 

(r + 1)^ ar+i(m, /i) = 2 + m(m — 1) r — yU ar(m, /i) — (r — 1)^ ar-i(m, /i) , 
(r + 1)^ br+i{m, /i) = 2 + m(m — 1) r — /i fe^l"^, lA ~ ~ 1)^ ^r-i(^5 /^) 



+ 



We choose, 



6 + m(m — 1) ar(m, /i) — 3 (r + 1)^ ar+i(m, yu) — 3 (r — 1)^ aj,_i(m, yu) . (102) 



ao(m, /i) = 1 , &o = , (103) 

and all the coefficients ar(m, /x) and hr{m, /i) become uniquely determined by eqs.( |102[ ). In 
particular, 



ai(m, jj) = —jj, , 8 a2{m, /i) = — /i [2 + m(m — 1) — /i] , 
6i(m, /i) = m(m — 1) + 3/^ , 

8 b2{m, /i) = [2 + m(m — 1) — /i] 6i(m, yu) — /i [6 + m(m — 1)] — 12 a2{m, fi) . 



(104) 



The coefficients of the leading pole singularities ai{m, fi) of /i and /2 obey identical recurrence 
relations as a consequence of the Baxter equation. 

Thanks to eq. (|102|) , (x; m, yu) for Q = /i 2 is an entire function of x. We find from eq. (l99|) 
that it has generically the form 

(x; m, /i) = [6 — m(m — 1)] lim [x Q (x; m, /x)] . 
The Baxter equation is fulfilled provided the limit in the r. h. s. is zero. 

Notice that a solution of the Baxter equation multiplied by a periodic function of x with 
its period equal to 1 is again a solution of the Baxter equation. An example of such a function 
is vrcotvrx which is a constant at infinity. Furthermore, if Q (x; m, yu) is a solution of eq. (pOj) , 
then Q (— x; m, — yu) is also a solution of eq.(p9|). 

We can now form linear combinations of the auxiliary functions /i (x; m, yu) and /2 (x; m, yu) 
in order to obtain true solutions of the Baxter equation [p!9[| : 



Q(^) (x; m, y(i) = /2 (x; m, yu) + B{m, fx) fi (x; m, yu) 



(105) 



Q(^) (x; m, y(i) = /i (x; m, /i) + C(m,y(i) /i (-x; m, -y(i) 



(106) 



As noticed above, the Baxter equation {x; m,fi) = is fulfilled at infinity provided the 
coefficients of in (x; m, /i) and in Q^^' {x; m, fi) vanish for large x, 

oo oo 

br{m, jj) + B{m, jj) ^ a.r(m, /i) = , 

r=l r=0 

oo oo 

aj,(m, /i) — C(m, /i) ^ ar(m, — /i) = . (107) 

r=0 r=0 

This gives for the coefficients the result, 

Therefore, the solutions Q^^-* (x; m, /x) and Q^^-* (x; m, /i) are completely determined. 

The functions Q^"^^ (x; m, /i) , Q*^^) (x; m,/i) and Q^^-* (x; m, /x) = Q*-^-* (— x; m, — /z) are re- 
lated by the linear equation 



[tt cot TTX + 5{m, /i)] Q*-^-* (x; m, /i) = Q''^'* (x; m, /i) — C(m, /i) Q*-^-* (— x; m, — /i) (109) 

where the coefficients in front of Q^^-* (x; m, /i) and Q*-^-* (— x; m, — /i) are calculated from the 
leading asymptotics at x — 1 and x —1, respectively. This is a special case of eq.(0) for 
n = 3 and r = 1. 

By finding the residues at single poles at x = 1 in the both sides of the above equation one 
can calculate 5(m, /i) 

C{m,fi)^ar{m,-fi) l^a,.(m,/i) m(m - 1) - 1 

5(m, /i) = B{m, fi) 2^ — 2^ 3 



fi r + 1 /i ^=2 ^ - 1 

The total energy for the Odderon described by the function Qm,m,fi being a bilinear combination 
of holomorphic and anti-holomorphic Baxter functions Q*-"-*, Q'^^^ and Q^^-* is given in ref.f 

d d 



E = i lim -— — — In 

A,A*^* dX dX* 



(A-z)2(A*-z)2|ArQ^,~,^(t)] , (110) 



Thus, the energy is expressed in terms of the ratio of residues for the single and double poles 
of the Baxter functions at x = 1. 

The holomorphic energies for the Baxter functions Q'-^-* (x, m, /i) and tt cot(7rx) Q^^^ (x, m, /x) 
should be the same, which leads to the quantization of fi in accordance with ref. [Il9| 

(5(m,/i)=0 (111) 

This equation fixes the possible values of fi for given m. 

Because in the bilinear combination m ^ ( ) there is a term proportional to the product 
of Q^"^^ (x; m,/i) and Q^"^^ {—x*; m, — /i), we obtain for the total energy 

E{m,m, fi) = B{m, fi) + B{Tn, —fi) . (112) 

We found for m = m = 1/2 the first roots numerically from the above equations (cf. |T9[) 



/ii = 0.205257506 .. . , /zg = 2.3439211 .. . , /xg = 8.32635 .. . , /i4 = 20.080497 .. . 

(113) 



with the corresponding energies 



El = 0.49434 . 



E2 = 5.16930. 



= 7.70234. 



Ea = 9.46283 , 



(114) 



We have followed the eigenvalues Ei, /xi as functions of m for < m < |. The result is 
plotted in fig. 0. Notice that only m = 0, 1 and | are physical values. For other m the curve 
describes the behaviour of the anomalous dimension for the corresponding high-twist operator 
(see sect. 5). 



0.5 



energy vs. m 
muvs. m 




Figure 1: The energy and fi as functions of m for the odderon eigenvalue Ei, /ii in the interval 
< m < 1. The picture is symmetric under m ■<=^ 1 — m. 

The energy vanishes at m = 0. This could be inferred from the fact that E{m = 0, /j, = 



0) = in the expression for E{m, = 0) given by eq.(74) of ref. ||T9| 



E{m, fi = 0) 



TT 



sm vrm 



+ ilj{m) + - m) - 2V^(1) . 



It should be noticed that E{m, ^ = 0) describes an eigenvalue providing that the function Q^^^ 
does not enter in the bilinear 
case 6{m, jj) does not vanish. 



does not enter in the bilinear combination of the total wave function ~ ,, because in this 



For the first eigenvalue we obtain numerically, 

Ei{m) "'=° 2.152 . . . m-2.754 . . . m^ + 0{m^) , /ii(m) ™=° 0.375 . . . v^-0.0228 m+0(mi) 

(115) 

Notice that all quantities are functions of m through the combination m(l — m). Therefore, 
they are invariant under the exchange m 1 — m. Thus, eqs.( |115| ) yield also the behaviour of 
El and /ii near m = 1 



m— »1 



2.152... d 



2.754... (1 -mf + O 



[m — 11 



/ii(mj 



m— >1 



0.375 ... VT^m - 0.0228 (1 



O 



m 2 



(116) 



The state with m = 1 and m = (or viceversa) is therefore the ground state of the 
Odderon. It has a vanishing energy and is situated below eigenstates (|114|) with m = m = 1/2. 

One can invert eq.( |116| ) to obtain the anomalous dimension 7 [see sec. 5] of the operator 
with m = I — 7 for ^ 0. The anomalous dimension is not small at arbitrary a/uj. On the 
other hand, we can not interpret the state with m = 3/2 as a physical state with = 2, = 
and a negative energy (see fig. 3), because is pure imaginary for this state according to fig. 4. 
Note, that next-to- leading corrections to the BFKL pair kernel (cf. |^) can give contributions to 
the eigenvalue u of the order of / (m — |), which will lead to the usual perturbative expansion 
of the anomalous dimension 7 ~ /uo for the corresponding operator. 

Furthermore, we consider the states with conformal spin n = 1 where 

m = 1 + iu , m = iv 

For small v we compute the energy of such state from eqs. (|115D and (|116|) with the result 

E = Eiim) + Eiim) = 5.51 ...u^ + O (u^) . 

The zero energy state that we find for u = is the one found in ref. |jl6| by a different approach. 
We have also followed the first eigenstate with n = as a function of u for real u and 



m 



1 
2 



The result is plotted in fig. ||. 
The eigenvalues for small v are 



Ei{v) = 0.49434 . . . + 1.8179 ... z/^ + 0(z/^ 



/ii(z/) = 0.205257 . . . + 0.48579 ... z/^ + 0{u^) 



We continue the first eigenstate with n = for m > 1 turning /i to purely imaginary. In 
figs. ^ and ^ the energy and Im fi as functions of m are plotted for the first eigenvalue in the 
interval 1 < m < 2. This describes the dependence of the anomalous dimension as a function 
of ac/uj for the corresponding operator with m = 2 — 7 (see the previous section). 



Near m = 2 the energy diverges while tends to zero. We find its behaviour near m = 2 
as follows 



EUm) — ^ + l + 2-m + 0\{m-2f 
Thanks to the m ^ 1 — m symmetry we have for m 



z/i 2 - m - - (m - 2)2 + C 
2 



:m-2) 



(117) 



^1 



m— > — 1 



m + 1 



+m+2+0 



[m + 11 



Now in order to compute the energy of the state with conformal spin n = 3 we use z/ — > as a 
regulator. We have from eqs. (|83D 

m = 2 + iv , m = —1 + iv . 

Hence, 

E = Ei{m) + Ei{m) 

Inverting eq.( |117] ) as discussed in sec. 5 yields for the anomalous dimension for 7 = 2 — m 



— + + 1 + W + Oip"^) = 2 + 0(iy^] 

IV IV 



0, 



7 



TX UJ 



- 2 



aNr 



TT UJ 



aNr 



TT UJ 



+ 



aNr 



TX UJ 



8 The solutions of the Baxter equation for four reggeons: 
the quarteton 

The Baxter equation for the quarteton (four reggeons state) takes the form 

i?4 (x; m, /i, q^) = 2a;^ + m(m — l)x^ — fix + Q {x; m, fi, q^) 

-{x + 1)^ Q {x + 1; m, /i, 54) - (x - 1)'^ Q {x - 1; m, /i, ^4) = . (118) 

where q^ = ifi , Im (/x) = 0. A new integral of motion ^4 appears here, and q^ are assumed to 
be real, which is compatible with the single-valuedness of the wave function in the coordinate 
space. 

Notice that if Q (x; m, /x, q^ is a solution of eq. ([L18|) then Q (— x; m, — /i, q^ is also a solution 



of eq.QTTS 



Following the general method presented in sec. Ill and ||T^ we seek solutions of the Baxter 
equation for the quarteton as a series of poles. We start by finding recurrence relations for the 
residues of the poles. Then, we impose the validity of the Baxter equation at infinity which 
gives further linear constraints on the pole residues. 

The auxiliary functions (^) take here the form 

ai{m, fi, qi) bi{m, n,q^) Ci{m, /i, q^) 



h [x] m, /i, 54) = 



=0 L 



(X - If 



(x - ly 



X — I 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 



Figure 3: The odderon energy as a function of m for the eigenvalue 1 in the interval 1 < m < 2 
and imaginary /i. 



i mu vs. m 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 

Figure 4: Im/i as a function of m for the odderon eigenvalue 1 in the interval 1 < m < 2. [Here 
Re/x = 0]. 



1=0 L 

00 

/i (a;; m, /i, ^4) = J2 



ai{m, n,q4) ^ bi{m, /i, ^4) 



1=0 



{x - If 
ai{m, fL,q4) 
X — I 



X — I 



(119) 



Imposing the Baxter equations at the poles x = I ,1 = 0,1,2,3,... yields the recurrence 
relations: 



(/ + l)^ai+i(m, /i,g4) 
(/ + l)^6;+i(m, /i,g4) 



2Z^ + m(m - 1) - yu/ + a;(m, yU, ^4) - (/ - 1)'^ a;_i(m, /i, ^4) , 
21^ + m{m — 1) f — fxl + q^ bi{m, fi, q^) — (/ — I)'* 6;_i(m, fx, q^) 



+ 



81^ + 2m{m — 1)1 — n ai{m, /i, q^) — 4 (Z + 1)^ aj+i(m, /i, 54) — 4 (/ — 1)'^ a;_i(m, /x. 



(/ + l)^Q+i(m, /i, ^4) = 2/^ + m(m - 1) - /i/ + 54 Q(m, /i, ^4) - (/ - l)^Q_i(m, /x, ^4) 



+ 2m(m - 1)/ -yu 6z(m, /x, ^4) - 4 (/ + l)^6;+i(m, /i, ^4) - 4 (/ - if bi^i{m, ^i^q^) 



-6(/ + 1)^ ai+i{m, yU, 54) - 6(/ - 1)^ aj-i(m, /i, 54) + [12/^ + m{m - l)]a;(m, /i, ^4) 



(120) 



We choose, 

ao(m, /i, 94) = 1 , 60 = , Co = , (121) 

and all the coefficients ai{m, fj,,q4,), bii^m, /x, 54) and Q(m, /x, 54) become uniquely determined 
by eqs.( p.20| ). In particular. 



ai(m, /i, ^4) = g4 , 61 (m, /i, ^4) = -4^4 - n 
Ci (m, /i, ^4) = 10 g4 + 4/i + m(m — 1) . 



(122) 



Taking linear combinations of fi (±x; m, ±/i, ^4) , /a (x; m, /i, ^4) and /s (x; m, yU, ^4) as in 
eqs.(^) and (^Dp we form three independent solutions Q^^^ of the Baxter equation ( [1181 ): 

Q*^^) (x; m, fi, qi) = /s (x; m, yU, ^4) + ai(m, yU, ^4) /2 (x; m, yu, 54) + 02 ("^, /i, ^4) /i (x; m, yU, ^4) , 
Q*^^^ (x; m, yu, 54) = /a (x; m, yU, ^4) + 7i(m, yu, 54) /i (x; m, yU, 54) + 

+ 72(m,y(i,g4) /i (-x; m, -yu,g4) . (123) 

The solution Q*^^) (x; m, yU, ^4) is proportional to Q*-^-* (— x; m, — yU, ^4). 

The recurrence relations for the coefficients of the poles guarantee that B4 (x; m, fi, q^) is 
an entire function. More precisely, we find from eqs. (|118|) and ( |119| ) that -B4 (x; m, yu, ^4) = 
ki + k2X, where ki and /c2 are some constants. These constants vanish and the Baxter equation 
is fulfilled provided the coefficients of x~^ and x~^ vanish for large x in the Q's of eq. (|123|) . 
Because the coefficients of x~^ and x~^ according to eqs.( p.31| ) from Appendix A are nonzero 
for the auxiliary functions fi (x; m, fi, q^) , /a (x; m, fi, q^) and /s (x; m, fi, ^4), these functions 
are not solutions of the Baxter equation. 

The coefficients ai(m, fi, q^), a2{m, fi, q^), 71 (m, fi, q^) and 72 ("^, /i, q^) are chosen imposing 
the Baxter equation at infinity. We present the linear equations on ai{m, fi,q4), a2{m, fi,q4) , 
■ji{m, fi,q4) and 72(771, yu, 54) and their exphcit solutions in Appendix A. 



There is one linear relation (0) among the Baxter solutions Q'-*-' for n = A 

S^'^^H, m, ^4) + TT cot TTx Q^"^^ {x; m, /x, ^4) = Q^^^ {x; m, /x, q4)+a^'^\n, m, g4)(5*-^'' (— x; m, — /x, ^4) 

where ^'■^''(/i, m, ^4) = m, ^4) — 7i(/i, m, ^4). A proof of this equation is given in Appendix 

B. Another relation is obtained from above one by the substitution 

The energy of the four Reggeons state is expressed through the bilinear combination m ii 



of the Baxter functions 



E = i lim — — In 



{\-iy\\* -iY\\\' Q 



m, m, /I, 54 



(^) 



dX dX* 

Again, E is related to the behavior of the Baxter functions Q^^^ near x = 1. That is, 



(124) 



.1 



TT cot[7rx] Q*-^-* (x; m, /x, 54) 



l-(x-l) U + + ai{m,fi,qM + O ([x - 1]' 



fx -11 



l-(x-l) |^4 + ^ + 7i(m,/i,g4)^ + O ([x - 1]= 



(125) 



We obtain for the total energy of the three solutions from eqs.( |124|) and ( |125|) , 

E(3)(m,m,/i,g4) = -ai(m, /i, ^4) - ai(m, -/x, ^4) , 

E^^^ (m, m, /i, ^4) = -71 (m, fi,q4) - 7i ("i, -/W, ^4) • 
The eigenvalue conditions 

5(^)(/i,m,g4) ~ ai(/i, m,g4) -7i(/i, m,g4) = , S^'^\-iJ, rn, q^) = 



(126) 



(127) 



guarantee that the holomorphic energy is the same for these three independent Baxter solutions. 



.(3) 



(m, /i, qi) = e^^) (m, /i, q^j = e^^^ (m, /i, ^4) . 



Eqs.( 127 ) fix possible values of fi and ^4 for given m. 

We find from the above equations for m = m = 1/2 the first roots numerically as. 



/i = 



g4 = 0.1535892... 



E = -1.34832. . . 



fi = 0.73833 , 
/X = 0. . . , 



g4 = -0.3703 , 



E = 2.34105. . . 



-0.292782... , E = 2.756624. . . , 



/x = 1.4100... , 94 = 0.73852... , E = 3.3581 .. . 



^ = 0... , 54 = 1.79992... 



, E = 5.67117... , 



/i = 0. . . 



^4 



-2.1857905... 



E = 6.2819490. . . 



:i28) 



The first eigenstate was reported in ref. |2^ where it is incorrectly identified as the ground state 
for the quarteton [see nota added]. 



m 



The 



In fig.0 the eigenvalue equations ( |127| ) are plotted in the /i, g4-plane for m 
curves intersect at the eigenstates ( |128| ). 

We find for the first eigenvalues with m = 0, m = 1 corresponding to = — 1 in eg . (|83D . [ 
n = 1 gives the same state with m ^m]. 



fx 







?4 



0.12167. . . 



E 



-2.0799... 



/i = 0.51214... 



-0.33288... 



E = 2.2007. . . 







-0.2905426 . . . 



E = 2.441210. . . 



(129) 



Therefore, the ground state of the quarteton with \n\ = 1 corresponds to m = 0, m = 1. 
Its energy, E = —2.0799 ... is below the lowest energy E = —1.34832 ... of the states with 
m = m = |. 

We plot in fig. ^ the eigenvalue equations ( |127| ) in the fi, g4-plane for m = 0, m = 1. The 
curves intersect at the eigenstates (|129|) . It should be noticed that the curves in figs. ^ and || 
turn to be qualitatively similar. The eigenstates with m = 0, m = 1 follow from those with 
m = m = I by an analytic continuation in m, m. 

We have followed the first eigenvalue as a function of m for < m < | . The result is plotted 
in fig. p. Contrary to the Odderon case, the energy eigenvalue does not vanish for m = 0. The 
energy decreases with m for < m < | and takes the value E = —2.0799 ... at m = while 
takes there the value 0.12167. . .. Further, fi vanishes for this eigenvalue for all m [see fig. ^. 

For m = 2 the lowest energy state of four reggeons goes to minus infinity and vanishes 
[gs is zero for all m in this state]. We find the following behaviour. 



El 



+2+2-m+C 



m 



[m - 2)^ 



94 



(m-2Y + 



[m 



2f 



(130) 



Thanks to the m ^ 1 — m symmetry we have for m 

4 



El 



m— 1 



m + 1 



+ 3 + m + (m+1) 



Now in order to compute the energy of the state with conformal spin n = 3 we use u 
regulator as in the odderon case. We find in this way. 



as a 



E = Ei{m) + Ei{m) ''=° 4 + 0{v') 
Inverting eq.( p.30| ) as discussed in sec. 5 yields for the anomalous dimension for 7 = 2 



■m 



7 = 4 

TT UJ 



+ 



aNr_ 



TT UJ 



The state with m = 3/2 (with n = 2, u = 0) can be considered as a physical ground 
state for the quarteton because for it the eigenvalue of q^ is real according to fig. 5. It has a 
large negative energy comparable with the energy of the BFKL pomeron constructed from two 
reggeized gluons. But to prove this conclusion it is needed to construct a bilinear combination of 
the corresponding Baxter functions to verify the normalizability of the corresponding solution. 




Figure 5: The as a function of m for the quarteton eigenvalue 1 in the interval < m < 2. 
We have for this eigenvalue = for all m. 
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Figure 6: The energy as a function of m for the quarteton eigenvalue 1 in the interval < m < 2. 
We have for this eigenvalue fi — ior all m. 
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Figure 8: The quarteton eigenvalue equations q;i(/x, 771,^4) — 71 (a*, 54) = and 
ai{—ii, m, Qi) — 7i(— m, Qi) — ior m — 0, m — 1 in the /i, g4-plane 



conformal spin 


conformal weights 


energy 

eigenvalues 

pomeron 


energy 

eigenvalues 

odderon 


energy 

eigenvalues 

quarteton 


n = 


m = m = 1/2 


-8 ln2 = -5.545177 


0.49434 


-1.34832 


n = 1 


m = 1, m = 








-2.0799 


n = 2 


m = 3/2, m = -1/2 


8(1 - ln2) = 2.45482 


imaginary fx 


-5.863 


n = 3 


m = 2, m = —1 


4 


2 


4 



TABLE 1. Conformal spins, conformal weights and corresponding lowest energy eigenval- 
ues for the pomeron, odderon and quarteton states. The odderon state with imaginary fi is 
discarded as non-physical. The pomeron and quarteton states with n = 3 have the same energy 
and are related presumably by the duality symmetry. 

As displayed in Table 1 the intercept for the quarteton ground state possessing conformal 
spin n = 2 is larger than that for the BFKL Pomeron. This result is not very suprising 
because four reggeized gluons clustered into two pomerons have an even larger intercept. A 
large intercept for the state with the conformal spin 2 may lead to such unphysical results as 
negative cross sections. But it is known that the unitarization of scattering amplitudes is not 
solved within a framework where the number of reggeons is fixed. 
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Nota Added: After completion of this work we have seen the preprint |^2[ studying similar 
problems. 

It is stated in refs. ||2l|, ^ that the ground state for three and four reggeons corresponds 

We show here that the lowest energy eigenvalue 
1 and n = 2, respectively. 



to vanishing conformal spin m — m = 0. 
corresponds in both cases to a non-zero conformal spin n 



It is claimed in ref. that the eigenvalues /z can be complex. In our opinion the reason 
for this disagreement with our results is, that the authors of ref. did not take into account, 
that in the Baxter-Sklyanin representation there are two different expressions for the energy in 
terms of Qm(A),Q~(A*) obtained at A ^ ±i. These two energy values coincide if and only if 
V'mm eigenfunction of the Hamiltonian. It leads to the constraint eq. (|68|) on the spectrum 
of the integrals of motion. 

We have computed the energies for the odderon states displayed in table 1 of ref.|j22|. States 
(0,2) and (6,2) coincide with our first two eigenstates 1 and 2 [see eqs.( |113D and ( |114| )] within 



the three figures precision of table 1 of ref.[^. States (2,2), (4,2), (8,2) and (10,2) possess 
complex fi and complex energy. Moreover, they fulfill the real part of our eigenvalue equation 
(|111|) up to O{10~^) while the imaginary part of 6{m, fi) turns to be of the order one. Our 
understanding is that these states with complex /i are not eigenstates of the hamiltonian. 

We furthermore computed the energies for the quarteton states displayed in table 2 of 
ref. [0]. States (2,0), (2,2), (4,0) and (3,3) coincide with our quarteton eigenvalues [see 
eq. (|128 )j 1, 3, 5 and 6, respectively. (Again, within the three figures precision of table 2 of 



ref. [^]). States (3, 1) and (4, 2) possess complex and complex energy. The real parts of E^^^ 
and E^"^^ computed from eqs.( |126| ) coincide up to O{10~^) while the imaginary parts of E^^^ 



and E^"^^ differ by amounts of the order one. Our understanding is again that these states with 
complex g4 are not eigenstates of the hamiltonian. 

A Asymptotic constraints on the solutions of the Baxter 
equation. The quarteton case. 

We impose here the Baxter equation at infinity on the solutions for four reggeons. As derived 
in sec. VI, the coefficients of and x~'^ in Q^^^ {x; m, /i, q^) and Q*-^^ {x; m, fi, q^) for a; oo 
must vanish. 

For large x the series ( |119| ) yield the following formal expansions, 

/3 (x; m, /i, q^) '''=°° — ^ Cr{m, /x, q^) + — ^ ^ [6,.(m, fi, q^) - r Cr{m, /i, ^4)] + O (ir) , 



r=0 y-^' r=0 



00 00 ^ -j^ 



/2 (x; m, /i, 54) = — X!^r-("^, /^,g4) + 7-^5I["^("^'Ai,?4)-'^&r("^, /U,g4)] + C> , 



za;;ro 1^^;;^ ^ \x^ 

fi (x; m, fi,q4) '''=°° —^0^(771, /i, ^4) - — -^^r a^(m, fi,q4.) + O ( — ] . (131) 

For Q^^^ (x; m, /x, q^) we find the conditions, 

00 00 00 

ai(m,/i,g4) "^Kim, /i, ^4) + a2("^, At, '?4) "^cirim, /i, ^4) = - ^ Cr(m, /i, ^4) , 

r=0 r=0 r=0 

00 00 

ai(m,/i,g4) XI K("^' ^.q^) - rhr{m, qi)] - a2{m, ^, q^) '^r ar{m, fi,q4) = 

r=0 r=0 

00 

= '^[rcr{m, fi,q4) - brim, fi,q4)] . (132) 

r=0 

Imposing the same constraint on Q*-^-* (x; m, fi, 54) yields, 

00 00 00 

7i(m,/i, 54) '^ar{m, /x, 94) - 72("^, A*, ^4) '^ar{m, -fi,q4) = - "^brirn, /i, 54) , 

r=0 r=0 r=0 

00 00 

7i(m,/i,g4) Xra,.(m, /i, ^4) + 72(m, /i, 54) ^ r a,.(m, -/i, ^4) = 

T-=0 r=0 

00 

= X K(m, /i,g4) - r6r(m, /i,g4)] (133) 



r=0 



Eqs.( |132| ) can be easily solved yielding, 



oo 



ai(m,/i,g4) = -<!^ra^(m, /i,g4)^c„(m, /i, ^4) 

ZAoi^m, /i, g4j [^^Q 

00 00 ^ 

- ^ a„(m, /i,g4) ^ [rc^(m, /i, ^4) - br{m, yU, ^4)] [ , 

n=0 r=0 ) 

(00 00 



a2{m, /i, ^4) = <! Y (^n{m, fi, q^) ^ [a,.(m, /i, ^4) - r ^^(m, /i, ^4) 



-^6r(m, /i,g4) a„(m, /i,g4) ^ (134) 

r=0 n=0 ) 



where 



Ac^(m,/i, ^4) = Y (in{m, fi,q4) ^ [rbr{m, /i, ^4) - a^(m, /U, ^4)]-^] ^r{rn, ^.q^) ^ n a„(m, /x, ^4) 

n=0 r=0 r=0 n,=0 

We analogously obtain from eqs. ( p.33| ) 



00 



ji(m,fi,q4) = <j Y (^n{m, -fi,q4) Y K("^,/i,g4) - r ^^(m, /x, ^4)] 

00 00 ^ 

- ^ 6„(m,/i,g4) ^ra^(m, -/i,g4) I , 

n=0 r=0 ) 

00 00 

72(m, /i, ^4) = — I Y (^n{m, /i, qi) ^ [ar(m, /i, q^) - r br{m, /i, ^4)] 

A^l^m, /i, g4j [^^Q 

00 00 

+ ^ 6„(m,/i,g4) ^ra^(m,/i,g4) I , (135) 

n=0 r=0 J 

where, 

00 00 00 00 

A^(m,/i,g4) = J2 an{m,iJ,,q4)J2rar{m, -/i, ^4) + X! "«("^' -/i, ^4) XI ^^("^""' ^4) • 

?i=0 r=0 n=0 r=0 

B Linear Relations among the solutions of the Baxter 
equation for the quartet on case. 

We present here a proof of the linear recurrence relations (0) for the quarteton case. The 
Baxter solutions Q^^^ (x; m, /z, ^4) and Q^"^^ {x; m, /x, ^4) are related by 

[«! (/i, m, ^4) — 7i (/X, m, ^4) + TT cot ttx] Q^'^^ {x; m, q^) 

= Q^^^ {x; m, /i, qi) - 72(/i, m, qA)Q^'^^ {-x; m, -/i, ^4) . (136) 

This is eg .(15^) forn = 4, r = 2 and 5^^''(/i, m, q^) = ai(/x, m, (3'4)— 7i(/i, m, ^4), a'^^^(/i, m, ^4) = 
— 72(/i, m^q^). Notice that Q*-^-* (x; m, /i, ^4) is proportional to Q^'^^ (— x; m, — /x, ^4) according 
to eq. (p3D and therefore eq. (|136|) is the only independent three-terms linear relation between 
Baxter solutions. 



In order to prove this relations we compute its triple, double and simple poles at x = I E Z. 
In the course of these calculations we use the fact, that vrcotTrx as a function of x has unit 
residues in all these poles. 

Moreover, we have from eqs. ( |119|) and (|123|) that 



(2) / . N x^_r>o ar(m, /i, ^4) br{m, /i, Qi) + 7i(/i, ; m, q^) ar{m, yU, q^) 



{x — ry 



X — r 



+Ar(m, /i, q^j + O {x - r) 

' {x; m, qi) = -72(/i, m, ^4) : \- Ar{m, fi, q^) + O {x + r) 

X + r 



(137) 



Inserting eqs.(|137|) into the Baxter equation (|118|) yields the recurrence relations 



(r + l)"^ Ar+i{m, fi,q4) = 2r'^ + m(m - l)r^ - yur + q^ Ar{m, ^4) - (r - 1)"^ Ar-i{m, /i, ^4) 



+ 



) + 2m{m — l)r — fi 



and 



brim, /i,g4) - 4(r + lfbr+i{m, fi,qi) - 4 (r - l)H,._i(m, /i, ^4) 
— 6(r + 1)^ ar+i(m, /i, 54) — 6(r — 1)^ ar-i{m, fx, ^4) + [12r^ + m(m — l)]ar(m, /i, ^4) 
+7i(/x, ; m, 54) I 8 + 2m(m - 1) r - /i ar(m, /x, ^4) 

-4(r + l)3a,+i(m, /i,g4)-4(r-l)3a^_i(m, /i,g4)} . (138) 

(r + 1)'^ A,,+i(m, /i, ^4) = 2r^ + m{m - 1) + /xr + ^4 A,,(m, yU, ^4) - (r - 1)"^ A,._i(m, fi, ^4) 
+72(/i, ; m, q^) | 8 + 2m(m - 1) r + /i ar(m, -/i, ^4) 

-4 (r + 1)3 (m, -fi, ^4) - 4 (r - 1)^ a^_i (m, -/i, ^4) } . (139) 

These recurrence relations can be solved in terms of the coefficients a^, br and using eq. ( |120| ). 
One has to take into account that Ar{m, fj,, ^4) and Aj.(m, /i, q^) do not obey the initial condi- 
tions dm]). We find, 

Ar{m, n, q^) = Crim, /i, q^) + 71 (/i, ; m, q^) br{m, fx, q^) + Ao(m, /i, ^4) a^(m, /x, ^4) 

i,.(m, /i, ^4) = 72(/^, ; 94) br{m, -//, ^4) + Ao(m, /i, ^4) ar(m, -//, ^4) (140) 

Using eqs. (|119|) , ( |123| ) and ( |14(]| ) we find that eg . (11361) holds at all of its poles. Therefore, 
the r. h. s. and the 1. h. s. of eq.( p.36| ) can only differ on an entire function. Taking into 
account the asymptotic behaviour of the Baxter functions we conclude that this entire function 
is identically zero. 



A r\ 
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